Nonstationary regime for quasinormal modes of the charged Vaidya metric 
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We consider non-stationary spherically symmetric n-dimensional charged black holes with varying 
mass m{y) and/or electric charge q{v), described by generic charged Vaidya metrics with cosmo- 
logical constant A in double null coordinates, and perform a comprehensive numerical analysis of 
the fundamental quasinormal modes (QNM) for minimally coupled scalar fields. We show that the 
"instantaneous" quasinormal frequencies exhibit the same sort of non-stationary behavior reported 
previously for the four-dimensional uncharged case with A = 0. Such property seems to be very 
robust, independent of the spacetime dimension and of the metric parameters, provided they be 
consistent with the existence of an event horizon. The study of time dependent Reissner-Nordstrom 
black holes allows us to go a step further and quantify the deviation of the stationary regime for 
QNM with respect to charge variations as well. We also look for signatures in the quasinormal 
frequencies from the creation of a Reissner-Nordstrom naked spacetime singularity. Even though 
one should expect the breakdown of our approach in the presence of naked singularities, we show 
that it is possible, in principle, to obtain some information about the naked singularity from the 
QNM frequencies, in agreement with the previous results of Ishibashi and Hosoya showing that it 
would be indeed possible to have regular scattering from naked singularities. 

PACS numbers: 04.30.Nk, 04.40.Nr, 04.70.Bw 
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I. INTRODUCTION 

The Quasinormal Modes (QNM) analysis is a central 
tool in the investigation of gravitational perturbations 
of stars and black holes. (For recent comprehensive re- 
views of the vast literature on the subject, see [l[ and [2[.) 
The QNM analysis of black holes is particularly relevant 
due to the possibility of observation of the ringdown sig- 
nals in gravitational wave detectors such as LISA, see 
for instance. Besides their relevance for future ob- 
servations, the QNM can also reveal important infor- 
mation about the structure and behavior of solutions of 
the Einstein equations. Many generalized solutions have 
been investigated, including non-asymptotically flat 0" 
0], charged (see for a comprehensive analysis), and 
time-dependent ones The AdS/CFT conjecture, 

in particular, has motivated many QNM analyses of gen- 
eralized black hole solutions According to the con- 
jecture, the QNM frequencies carry information about 
the thermal properties of the associated conformal field 
theory 

The Vaidya metric [l^ was originally proposed to de- 
scribe the spacetime outside a radiating star. It has also 
been the usual starting point for the study of time de- 
pendent black holes QNM [9l-[l3j . Namely, it corresponds 
to a time dependent solution of Einstein equations for a 
spherically symmetric body immersed in a unidirectional 
radial null-fluid flow. It has been also widely used in the 
analysis of spherically symmetric collapse and the forma- 
tion of naked singularities for many years. (For further 



references, see [T^l and the extensive list of [1^). It is also 
known that the Vaidya metric can be obtained from the 
Tolman metric by taking appropriate limits in the self- 
similar case[l^. This result has shed some light on the 
nature of the so-called shell-focusing singularities [l^l , as 
discussed in detail in [l8l - l23| . The Vaidya metric has 
also proved to be useful in the study of Haw king radia- 
tion and the process of black- hole evaporation [24| - l27| , in 



28{ , and in recent numer- 
The charged version of 
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the stochastic gravity program 
ical relativity investigations [2S 
the metric is also an usual starting point to the study 
of many aspects of charged black hole physics and naked 
singularities j30l - l33 |. 

The main purpose of the present paper is to consider 
the QNM of time dependent backgrounds corresponding 
to a general Vaidya metric, with special emphasis on non- 
stationary effects, generalizing and further developing the 
work done in [ll|, [13] . This kind of QNM analysis is cer- 
tainly relevant from the physical point of view since, for 
instance, it is quite natural to expect that a real black 
hole be affected by processes which can change its mass 
such as, for instance, mass accretion or even Hawking 
radiation, which would indeed imply a decreasing mass. 
Any signal coming from a black hole could, in principle, 
have some nonstationary component. 

As a secondary objective, we show that our numeri- 
cal setup can be used to investigate the QNM frequen- 
cies in the case of a Vaidya metric evolving towards a 
Reisser- Nordstrom naked singularity. One should expect 
the breakdown of our approach in the presence of a naked 
singularity, of course. However, we show that it is in prin- 
ciple possible to obtain some information about the sin- 
gularity from the QNM frequencies, in agreement with 
the results of Ishibashi and Hosoya [35|, who conclude 
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that it is indeed possible to have regular scattering from 
naked singularities. 

This paper has three more sections. In the next one we 
present the main equations describing the most generic 
Vaidya metric in double-null coordinates and write the 
scalar wave equation in a canonical hyperbolic form, spe- 
cially suitable for the numerical analysis presented in Sec- 
tion III. In the last section, we present some concluding 
remarks. 



where ka is a radial null vector and h{u, v) is a smooth 
function characterizing the null-fluid radial flow. We will 
consider here, without loss of generality, the case of a 
flow along the w-direction. 

From the Einstein-Maxwell equations with metric ([2]) 
and energy-momentum tensor ^ wc obtain the following 
equations for the functions /, h and r: 



(4) 



II. THE SCALAR WAVE EQUATION IN THE 
VAIDYA SPACETIME 



The n-dimensional Vaidya metric was first discussed 
in (36| . It can be easily cast in n-dimcnsional radiation 
coordinates (w, r, 0i, . . . , 0„_2) as done, for instance, in 
■ The n-dimensional charged Vaidya metric in radia- 
tion coordinates, obtained originally in [sif . reads 
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where n > 3, c ~ ±1, and dD,^_2 stands for the metric 
of the unit {n ~ 2)-dimcnsional sphere, assumed here to 
be spanned by the angular coordinates O2, ■ ■ ■ , dn-2) 
in the usual way. For the case of an ingoing radial flow, 
c = 1 and m(v) is a monotonically increasing mass func- 
tion in the advanced time v, while c = — 1 corresponds 
to an outgoing radial flow, with m{v) being in this case a 
monotonically decreasing mass function in the retarded 
time V. The constant q corresponds to the total electric 
charge. In principle, one can also consider time depen- 
dent charges q as done, for instance, in [H]. This situ- 
ation will of course require the presence of charged null 
fluids and currents, whose realistic nature we do not ad- 
dress here. We have already reported some preliminary 
results on QNM for this case in [H] . 

We will now extend the approach proposed in [s^, li^l 
and derive the double-null formulation for the most gen- 
eral Vaidya metric: n-dimcnsional, in the presence of a 
cosmological constant, and with varying electric charge. 
Only the main results are presented. The reader can get 
more details on the employed semi-analytical approach in 
[39I lioj and the references cited therein. We recall that 
the n-dimensional spherically symmetric line element in 
double-null coordinates {u,v,0i, . . . , dn-2) is given by 



ds^ = -2f{u,v)dudv + r^{u,v)dnl_ 



■2i 



(2) 



where f{u, v) and r(u, v) are non vanishing smooth func- 
tions. The energy- momentum tensor of a unidirectional 
radial null-fluid in the eikonal approximation in the pres- 
ence of an electromagnetic field Fab is given by 
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where B(v), m{v) and q{v) are arbitrary integration func- 
tions. If we choose B{v) = constant, we can interpret 
m{v) and q{v) as the mass and charge of the solution, re- 
spectively. These two functions must be monotonic, and 
must be chosen in a way that satisfies the null-energy 
condition [s^]. The details of the derivation of eqs.Q- 
(|n]) and the full analysis of the charged Vaidya metric in 
duble null coordinates will be presented in another paper 
currently under preparation. 

For our QNM analysis we will consider the evolution 
of a massless scalar field governed by the Klein- Gordon 
equation 



= 0, 



(7) 



in the background Taking advantage of the spheri- 
cal symmetry, we decompose the scalar field in terms of 
higher-dimensional spherical harmonics 

^if = J2^{u,v)Ye^iei,...,9n^2), (8) 



for which 



v^y,„ = -^(£ + n-3)r,„, (9) 



where stands for the Laplacian operator over the (n- 
2)-dimensional unit sphere. Substituting the ansatz ([S 
in ([7]), we obtain for the metric 1^ 
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Using now the substitution ip = r we get 

'P,uv + V{u,v)(p = 0, (11) 

where 



^(^ + n-3)^ (n-2)(n-4) (n - 2) 
V = / — r ,,r_„ r , 
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The wave equation is already written in a canon- 
ical hyperbolic form and it will be the starting point for 
our QNM analysis. Nevertheless, the potential ([T2l) can 
still be cast in a more convenient way. From equations 
dH) and ©, we find 
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Now we can finally write the potential V{u, v) as 
'e{i + ?i-3) (n-2)(n-4) 
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The conventions we adopted for the mass, electric charge 
and cosmological constant appearing in the potential ([15]) 
and used in the derivations of this section are the stan- 
dard ones employed in the Vaidya metric literature. In 
particular, our expressions coincide, in the appropriate 
limits, with the previous results for the A = [31| and for 
the q = [1^ cases. However, the commonly employed 
conventions in the QNM literature are slightly different, 
see, for instance, [41|. 



III. NUMERICAL RESULTS 

The wave equation (jlip is already written in a canoni- 
cal hyperbolic form and it can be numerically integrated 
by means of a characteristic scheme, i.e., one can evolve 
Lp(u,v) in V knowing (p{u,vq). In order to determine 
V{u,v), one needs to known f{u,v) and r{u,v) and, 
consequently, Eq. ^ must be solved in each evolu- 
tion step. Since only ip{u, v) is required to determine 
(p{u, V + dv), one can implement the algorithm in an effi- 
cient way avoiding unnecessary calculations. We use here 
the same numerical scheme used in [Til . [l2j to evolve the 
system and read the QNM frequencies, which allows us 
to perform an exhaustive numerical analysis with modest 
computational resources. 

All mass functions considered in this work are of the 
form 

2m{v) = {mf + mi) + {ruf — rrii) tanh Pm{v — Vm), (16) 

where rrii and mf stand for the initial and final mass, re- 
spectively, and pm controls how fast the change is, with 
the maximum rate of change at u = Vm. We use an anal- 
ogous expression for the time dependent electric charge 
q{v). Our main results do not depend on the exact form 



of the functions m{v) and q{v). Choices such as (|16p . 
however, are very convenient due to their smooth behav- 
ior and the "asymptotically static" limits for v — > ±oo. 

For the numerical calculation of the u-evolution, the 
algorithm requires the evaluation of (p{u, v) for all values 
of u corresponding to the exterior region of the black hole, 
let us say, for Uh < u < oo, with u = corresponding to 
the event horizon. Technically, it is easier to control the 
approach to the regions near the horizon by introducing 
an appropriated kind of tortoise coordinate U such that 
the external region Uh < u < oo will correspond to — cx) < 
U <oo. 



A. Static case: Reissner-Nordstrom black hole 

It is instructive to start by testing our code with the 
standard Reissner-Nordstrom (RN) case, for which a vast 
set of results is available in the literature • The poten- 
tial (|15p for n — 4 and A = becomes simply 



V{u,v) 



2 V 



2 m 



2q2 



(17) 



In order to compare the QNM frequencies obtained from 
the numerical integration of (|lip for the potential p7)) 
with well-known results for the RN black hole, we per- 
form the following change of variables 



U = U-- 



2ri 



- — In ( r ( 

r_ V 2 ^ 



2r' 
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r_ V 2 



where 



(18) 



± 



(19) 



as usual for the RN black hole. This is necessary in order 
to guarantee that the QNM frequencies be defined with 
respect to the same scales. Moreover, the new tortoise 
coordinate U given by ([T5| allows a better description of 
the exterior region of the black hole, notably of the re- 
gions very close to the horizon r+. The numerical grid is, 
therefore, written inU xv coordinates, and u is obtained 
by inverting the definition psp . 

A comparison between the the numerical results ob- 
tained with our code and the values found in the litera- 
ture is shown in Fig. [T] We observe a very good agree- 
ment with the known values for the RN QNM frequen- 
cies. This simple example allows us to calibrate all the 
algorithm parameters in order to attain a pre-established 
accuracy. 



B. Time-dependent case 

We perform a comprehensive QNM analysis for the 
wave equation (|TT|) with the potential (|T5|) . We are 
mainly interested in the nonstationary regimes like the 
one described in [ll|, corresponding to the case q = 0, 
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FIG. 1. Frequencies of the quasinormal modes of a scalar 
perturbation with £ = 2 obtained for a RN black hole with 
our code (smooth line) compared with values found in the 
literature [ill,!!!]. 



n = 4, and A = in the potential ([T5|) . As we show 
in this section, the stationary and nonstationary regimes 
for QNM are also present for the general situation cor- 
responding to the potential ([15]). Moreover, we verify 
appreciable deviations from the stationary regime when- 
ever |m"(w)| or is larger than |a;/|, where w/ is 
the imaginary part of the first QNM frequency. This 
behavior occurs irrespective of the other parameters of 
the potential ((T5|) . provided they be consistent with the 
existence of a black hole. 

The case of an asymptotically flat (A = 0) "time de- 
pendent RN black hole", i.e., a RN black hole with time 
dependent mass m{v) and/or charge q{v) functions, is 
particularly interesting. The determination of the new 
tortoise variable U equivalent to (|18p is quite more in- 
volved for this case. The problem here is the definition, 
and the numerical determination, of the event horizon. 
(See [1^ for a discussion of the implications of this prob- 
lem in the semi-analytical approach used here.) 

For our purposes here, the event horizon r_|_ is the last 
null geodesic (up to the machine precision) escaping to- 
wards infinity, requiring the full numerical solution of (jH]) 
prior to the analysis of the wave equation pT|) . The de- 
termination of the real Cauchy horizon r_ is easier since 
any null geodesic inside the event horizon tends to r_ 
along the w-evolution. 

In Fig. [5J we present the time variation for the real 
and imaginary parts of the QNM of the scalar perturba- 
tions with £ = 2 for a 4-dimcnsional RN black hole with 
time dependent mass and charge functions. The QNM 
frequencies change in a similar way with the time vari- 
ation of the mass or the charge, but they are typically 
more sensitive to the charge variations. 

The results from the lower plot of Fig. [2] can also be 
seen in Fig. (jS]), in the LOr x i^i plane, where we see how 



the variation goes from one stationary state to another 
through a nonstationary trajectory. In the right plot of 
FiglSl we use the minimum value of the variation of uji 
in the nonstationary trajectories of the left plot to quan- 
tify the deviation from the stationary trajectory. These 
values are plotted against p, which can be related to a 
measure of how fast the variations are. We can see that 
the behavior can be described extremely well by a lin- 
ear fit at first, deviating from this fit for p > 0.1 (faster 
variations) . 

The case of asymptotically extremal (qf ~ nif) black 
holes is specially notable in our approach due to numeri- 
cal technicalities. Our results for this case are presented 
in Fig.m In the right plot we present a convergence test, 
essential in this case, with the results obtained for dif- 
ferent resolutions, while in the left plot we show how the 
results depend on the manner in which q{v) — > m. We 
remind here that an extremal RN black hole is, of course, 
yet a regular black hole and not a naked singularity, and 
its QNMs are well defined. It is only due to numerical 
limitations that most works in the literature do not reach 
this limit. 

The left plot of Fig. |4] shows once again that the devia- 
tion from the stationary behavior increases as we consider 
faster variations of the background, which are quantified 
by p. As this deviation becomes smaller, however, it is 
possible to note a nontrivial trajectory of the complex 
frequencies in the plane (more easily discernible in the 
right plot), whose origin is still unclear. 



C. Formation of a naked singularity 

We can use our formalism and numerical setup to 
probe the formation of a naked singularity in the space- 
time. The idea that it would be possible to have reg ular 
scattering from naked singularities is rather old [43 ■ In 
fact, as it was first noticed by Gibbons [3], minimally 
coupled scalar fields have the remarkable property of be- 
ing regular at the origin of a Reissner-Nordstrom solu- 
tion, where the spacetime manifold is irremediably sin- 
gular. The non-singular behavior of scalar [4^ and other 
[46| fields around a Reissner-Nordstrom naked singularity 
has been investigated since then. Despite being mathe- 
matically well posed, the field dynamics near a naked sin- 
gularity are typically ambiguous from the physical point 
of view since it is not clear which boundary condition one 
needs to impose for the field at the singularity. 

The main motivation of this part of our investigation 
is the work of Ishibashi and Hosoya [s^ , which demon- 
strated explicitly that is indeed possible to have unam- 
biguous and regular scattering from naked singularities. 
After all, the presence of a naked singularity might be 
less harmful than originally conceived. This has been ex- 
plored recently in a quantum field theory scenario in the 
series of papers [43-[50|. 

We apply our numerical setup to simulate a charged 
black hole losing mass while keeping its charge constant. 
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FIG. 2. Upper plot: LUr (left) and uji (right) as function of v for the £ — 2 scalar perturbation of a black hole with constant 
charge q = 0.25 and m{v) given by (|16[) with rrn — 0.5, m/ — 0.65, Vm = 75 and different values of pm- Lower plot: tur (left) 
and iJi (right) as function of v for the £ = 2 scalar perturbation of a black hole with constant mass m = 0.5 and q{v) given by 
(|16l) with Qi — 0.35, qf = 0.4, Vm = 75 and different values of pq. In all graphics, the horizontal lines show the frequency values 
for static RN black holes with the initial and final configurations. The fit for the frequencies is typically more sensitive for the 
values of oji. 
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FIG. 3. Left plot: ujr x LUi plane, showing the transition between equilibrium states for a varying charge case according to (|16|l . 
The values of pq for these curves are 0.04, 0.06, 0.08, 0.1 and 0.2 (from top to bottom). Right plot: the minima of oJi obtained 
in the transition shown in the left plot for different values of pq, to quantify the non-adiabatic behavior. 



leading to a Reissner- Nordstrom naked singularity. For a irif < q. At the point v^, when m{v,,) = q the evaluation 
fixed charge value q, we set a mass function ()16|) between of ip{u, u«) is problematic since for the case of a naked 
two values rrii and m f , and take the final mass to be singularity the full range of u must include the singularity 



6 




-0.11 I ' ' ' ' ' ' ' 1 -0.102 I ' ' ' ' — 

0.56 0.58 0.6 0.62 0.64 0.66 0.68 0.7 0.55 0.57 0.59 0.61 0.63 

mto f mo) f 



FIG. 4. Left plot: ujr x uji plane, showing the transition between equilibrium states for a varying charge case according to (|16|l . 
The final state is an extremal RN black hole (g = m) and the values of pq are the same of Fig. O Right plot: Convergence test 
made for a fixed value of pq = 0.04 and three different resolutions obtained with varying values for the integration stepsize h. 
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FIG. 5. Left plot: mass functions of the form (|16|) with rrii = 0.65, pm = 0.1, Vm = 75 and different values of rrij, used 
to force the formation of a naked singularity while keeping q = 0.5 constant. Right plot: ujr as function of v for the 1 = 2 
scalar perturbation of black holes with constant charge and varying masses as given in the left plot. The horizontal line shows 
the approximate frequency for an extremal RN black hole. The behavior of oj^ is quite similar but typically known with less 
precision. 



at r = 0. Despite the fact that the scalar field can indeed 
be finite at r = |45|, our code cannot deal properly 
with the terms like V {u,v)i^{u^v) in the vicinity of the 
singularity. This point is now under investigation, and a 
better numerical code is being developed for this case. 

Our objective here is to investigate how the QNM be- 
have until the very last moment before the formation 
of a naked singularity. The results obtained for the a;^ 
part of the QNMs are presented in Fig. ([S])- AH simula- 
tions end at the point when an extremal RN black hole is 
formed, but the frequencies are all different! This is the 
result of the different nonstationary trajectories followed 
by the QNM in each case. We remark here that in the 
cases with a more rapid variation m'(w*) it seems there 
is not enough time before the simulation stops to see the 
maximum deviation of the frequency and the subsequent 
relaxation to a possible asymptotic value. 



IV. FINAL REMARKS 

Our results confirm, for the general Vaidya case, the 
same nonstationary behavior corresponding to the inertia 
of the QNM frequencies identified in In particular, 
for situations with r" (v) > a;i(u), the QNM frequencies 
will not follow in the uj^ x uji plane a trajectory corre- 
sponding to the instantaneous frequency associated to 
a RN black hole of a given qj-m ratio (see Figs. [3] and 
2]), and the inertial behavior of the QNM can be iden- 
tified. Moreover, we see that the faster the change in 
m{v) or q{v) is, the bigger the deviation from the sta- 
tionary regime will be. We could determine that the 
deviation from the stationary behavior is proportional to 
p for masses and charges varying according to (1161) . at 
least for values of p < 0.1, see Fig. |31 

Similarly to the uncharged case (Tl| . the imaginary 
parts of the QNM frequencies are typically more sen- 
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sitive to mass and charge changes than the real parts. 
However, we have identified a new behavior whose physi- 
cal origin is still unclear for us: the QNM frequencies are 
typically more sensitive to charge than mass variations 
(compare the upper and lower graphics in Fig. [2]). 

This leads to a curious fact and a further conclusion. 
Consider two situations such that, in the first, electric 
charge is constant and mass is increasing, let us say ac- 
cording to ^TE\\ . In the second, mass is constant and 
electric charge is varying in a way that both cases have 
the same quotient q/m as function of w. For these two sit- 
uations, the instantaneous QNM frequencies will evolve 
differently in the cOr x uji. According to our results, the 
case of varying charge will typically show a bigger de- 
viation from the stationary curve q/m. This new be- 
havior demonstrates explicitly that the ratio q/m is not 
enough to characterize a non-stationary black hole. (This 
does not represent any challenge to the no-hair theorems, 
though, since they typically deal with stationary solu- 
tions.) 

As for the q — case considered in [ll| , an interesting 
extension of this work would be analysis of the highly 



damped QNM. Since for such overtones the ratio \uJi/ujr\ 
is always larger than for the fundamental (n = 0) QNM 
considered here, including, for sufficient large n, cases for 
which jwi/wrl > 1, it would be interesting to check if the 
stationary behavior could be improved for n > 0. In par- 
ticular, it would be very interesting to test the high sen- 
sitivity of the instantaneous QNM overtones to electric 
charge variations. We remark that the numerical analysis 
presented here cannot be directly extended to the n > 
case since one cannot identify the overtones numerically 
with sufficient accuracy. We believe, however, this could 
be attained by means of the WKB approximation. 
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